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STABILITY OF GAS MEASURES UNDER PERTURBATIONS AND 

DISCRETIZATIONS 


ROBERTO FERNANDEZ, PABLO GROISMAN AND SANTIAGO SAGLIETTI 

Abstract. For a general class of gas models —which includes discrete and continuous 
Gibbsian models as well as contour or polymer ensembles— we determine a diluteness 
condition that implies: (1) Uniqueness of the infinite-volume equilibrium measure; (2) 
stability of this measure under perturbations of parameters and discretization schemes, 
and (3) existence of a coupled perfect-simulation scheme for the infinite-volume measure 
together with its perturbations and discretizations. Some of these results have previously 
been obtained through methods based on cluster expansions. In contrast, our treatment 
is purely probabilistic and its diluteness condition is weaker than existing convergence 
conditions for cluster expansions. 


1. Introduction 

Phase transitions in statistical mechanics are often stndied throngh sequences of models 
involving convergent sequences of parameters. The analysis is usually “towards a target”: 
properties of the target model are inferred from properties of models in the sequence. 
For instance, sequences of models with asymptotically vanishing helds can be used to 
show hrst order phase transitions for the zero-held model. In the same spirit, transitions 
in continuum-space models are often studied through limits of models on lattices with 
decreasing mesh. 

In this paper we treat the opposite —“from the target”— point of view. We investigate 
conditions under which properties of a target model are inherited by models obtained 
by perturbing parameters (e.g. helds or fugacities) or the conhguration space itself (e.g. 
through discretizations). While our treatment is general, our basic motivation came from 
target models in the continuum, for which we wished to address two types of issues: 

Faithfulness of simulation schemes: Simulation of continuum models requires 
unavoidable discretizations. It is tacitly understood that using sufficiently rehned 
discretizations leads to trustable determinations of phase diagrams. Still, one may 
wonder if this is always the case. The question involves, in fact, a previous issue: 
How should the discretization be performed? The natural choice would be to 
discretize the model (both conhguration space and interactions) and sample from 
the Gibbs measure in the discrete system. An alternative, however, would be to 
sample the discretized version of the actual continuum measure. These approaches 
are quite likely not equivalent in general, as it is known that coarse-graining may 
lead to non-Gibbsianness. 

Universality: Stability of the continuum equilibrium measure under discretizations 
should imply the irrelevance in the discretized systems of interaction terms which 
disappear in the continuum limit. This includes, for instance, hard-core exclusions 
involving events with zero probability in the continuum (see the thin rods model 
in Section [2.3.51 below, for example). In this sense, stability can be interpreted as 
the continuum model acting as a “universality class” for whole families of discrete 
systems. 
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Our treatment is geared towards general gas models, that is systems involving families of 
geometrical objects —possibly with further decorations such as color or spin— distributed 
on some underlying space. This distribution is assumed to be, while on bounded volumes, 
absolutely continuous with respect to a basic “free” Poisson process. Our approach uses 
only general properties of this density with respect to the “free measure”, and hence it is 
applicable to general point processes not necessarily endowed with a Gibbsian description. 
In particular, it applies to the contour ensembles used to describe low-temperature phases, 
starting with the well-known Peierls contours. 

The results reported below hold for models satisfying an appropriate dilution condition 
which physically corresponds to gaseous phases. This dilution condition leads to objects 
typically clustering into hnite islands separated by percolating empty space. In particular, 
our condition implies uniqueness of the inhnite-volume point process or Gibbs measure. 
In the Gibbsian setting, our requirement leads naturally to models at high temperature 
or low fugacity. However, the generality of our approach makes the technique relevant 
also for low-temperature (or condensation) regimes in which typical conhgurations can be 
described by diluted contours. This generality will be exploited in a forthcoming paper; 
in the present paper we focus on stability issues within the uniqueness regime. 

Our main result establishes that whenever a target model satishes this dilution condition 
all sufficiently small perturbations of it also admit exactly one consistent inhnite-volume 
measure which, furthermore, is a slight perturbation of the measure in the target model. 
Let us list some particularities of our approach: 

• The dilution condition amounts to a strong form of uniqueness that is, however, 
weaker than the condition associated to the validity of cluster expansion methods 
and therefore applies to a wider range of systems. This extension comes at a cost: 
Roughly speaking, our condition implies only the continuity of the unique measure 
with respect to slight perturbations in the parameters of the model, whereas 
the convergence of cluster expansions leads to analytic dependences. 

• Our approach yields a coupled perfect simulation algorithm, that is an algorithm 
yielding simultaneously exact samples of the inhnite-volume equilibrium measures 
(restricted to a hnite volume) of target and approximating models. 

• This coupled algorithm leads to the almost sure convergence of the samples above, 
and thus to the weak convergence of the equilibrium measures of perturbed models 
towards that of the target model. 

We believe that our results are quite natural and easy to apply, as we illustrate through 
a number of examples. Nevertheless, our presentation is not devoid of technical details. 
In particular, in Section [2] we present a careful account of the general setup for gas models 
(which is used indistinctly for discrete and continuum systems) followed by the precise 
dehnition of “approximation” operations (i.e. perturbations of the conhguration space). 
Also, in Section [5] we discuss with some detail the ancestor algorithm which constitutes 
the main tool of our analysis. This algorithm —introduced 15 years ago as a substitute 
for cluster expansions in [S]— reconstructs conhgurations through a (time-backwards) 
oriented percolation model of space-time “cylinders”, i.e. objects in the gas model which 
live for a certain period of time. This ancestor algorithm succeeds —implying uniqueness 
and space-time mixing of the inhnite-volume measure— if these cylinders do not percolate. 
Such a condition is naturally suited for stability studies, because hnite cylinder clusters are 
robust under perturbations. Our dilution condition is crafted to ensure, in general terms, 
the lack of percolation for cylinders associated to the target model. This lack of percolation 
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is then inherited by the perturbed models, whose cylinder clusters become in one-to-one 
correspondence with those of the target model as the strength of the perturbation vanishes. 

2. The basic setup 

2.1. Configuration space. We start by describing the general measure-theoretical setup. 
The dehnitions aim at a general conhguration space on an underlying space of locations. 
Whenever the latter is discrete, we call the model in question a lattice system, although 
other type of systems also fit into this framework. 

2.1.1. Particle configurations and configuration space. We consider two locally compact 
complete separable metric spaces: a location space {S, ds) and a space of animals {G, do). 
A countable S is often called a lattice and a hnite G is interpreted as a set of colors or spins. 
The product space S x G is also locally compact complete and separable if endowed with 
the metric d = ds + dc- For convenience, we shall call an element (x, 7 ) G S' x G a particle 
and denote it simply by We interpret it as an animal 7 positioned at location x. 

The general dehnition of conhguration space requires special care. For our purposes, 
it will be convenient to adopt the general framework of point processes featured in m- 
In this framework, conhgurations are identihed with locally hnite measures on S' x G 
obtained as a superposition of delta-measures signaling the presence of particles. 

Notation 2.1. Given a metric space {X,d) we denote: 

• By Bx the class of all Borel subsets of (X, d). 

• By Bx the set of elements of Bx with compact closure. 

We recall that a set i? G Bx is locally finite if for every B' G Bx the set i? Hi?' is hnite. 
Also, a measure ^ on {X, Bx) is called a Radon measure if ^{B) < -|-cxo for every B G B\. 
Conhgurations correspond to particular Radon measures supported on locally hnite sets. 

Definition 2.2. A Radon measure f on (S' x G, Bsxg) is said to be a particle configuration 
if e(5) G No for every B G B^^g- 

The following proposition states that a particle conhgurations can be actually identihed 
with a locally hnite collection of particles, in which particles may appear more than once. 

Proposition 2.3 ([HI Lemma 2.1]). A measure on S' x G is a particle conhguration if 
and only if there exist a locally hnite set ^ S x G and a map : (.^) —)■ N such that 

( 2 . 1 ) ^ 

7x6(0 

with 6^^ the Dirac measure centered at 7 ^,. We call fl 2 .ip the standard representation of 

Thus, a particle conhguration ^ can be thought of as a locally hnite family of particles, 
each with its own respective multiplicity. The set (^) shall be called the support of 
It is the collection of particles which constitute with no regard for their multiplicities. 
If we take these into consideration then we obtain the weighted support of ^ given by 

[^] := {hx,i) G (S' X G) X N : 7a; G (0 and 1 < f < rn^(7x)}. 

Dehnition l2.2l alone is sufficient to dehne the conhguration space whenever G is compact. 
For the general case, however, we will require an extra restriction. 

Definition 2.4. 

i. A measure ^ on (S' xG, Bsxg) is said to be of S-locally finite allocation if it satishes 
^(A X G) < -foo for every A G B^g. 
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ii. The configuration space of S' x G is the space J\f{S x G) of all particle configurations 
on S' X G which are of S'-locally finite allocation. 

iii. The set of configurations (supported) on A G is the set 

AA(A X G) := e fif{S X G) : (0 C A X G}. 

[The notation is slightly abusive.] 

We now define the natural notions of restriction and superposition of configurations. 

Definition 2.5. Given ^ G A/”(S' x G) and A G BsxG, fhe restriction of ^ to A is defined 
as the particle configuration such that for every B G Bsxg 

UB) = aAnB). 

Equivalently, if f := then is given by the standard representation 

7:.G{0nA 

To improve readability, in the following we will write instead of ^axg for A G B^g. 

Definition 2.6. The superposition of two configurations cr, r; G J\f{S x G) is defined as 
the particle configuration a ■ rj such that for every B G Bsxg 

{a-v)iB) = aiB)+piB). 

In the particular cases in which a G J\f{A x G) and i] G x G) for a certain A G Bg, 

the superposition a-rj can be thought of as a concatenation. In these cases, the operations 
of restriction and superposition induce a natural identification between N'{S x G) and 
AA(A X G) X AA(A'^ X G) for any given A G B'g. Indeed, the applications 

U{S X G) ^ M{A X G) X A^(A" x G) 
k '-1 (^A,^A'=) 

and 

N{A X G) X AA(A" x G) ^ N{S x G) 

{(T,r]) I—^ a-7] 

are bijections and have each other as their respective inverse. 

2 .1.2. Measurable structure. The space J\f{S x G) is endowed with a measurable structure 
by considering the a-algebra B generated by the counting events, i.e. 

(2.2) jr = a ({^ G X G) : i{B) = A:} : A: G No and 5 G B%^g) . 

Alternatively, if one considers for every B G BgxG fhe counting variable 

Nb : M{S X G) ^ No 

AfaW = niB). 

then 

7-= ff (Affi : B e BLc) . 

More generally, for any A G BgxG fhe a-algebra Ba of events occurring in A is defined as 
the one generated by the counting events inside A, i.e. 

Ba := a{{f eAf{SxG):f{B) = k} :ke No and B eB°^) 

= a{NB:BeB^^). 

The case A = A x G for A G is of particular relevance. First, we make the following 
important observation. 
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Remark 2.7. The identification between J\f{S x G) and J\f{A x G) x x G) defined 
above is in fact a measurable isomorphism if the spaces are endowed with the a-algebras 
and J^axg J^a<=xg, respectively. 

Further, we introduce the natural notions of local events and observables. 

Definition 2.8. 

i. A function / : J\f{S x G) —)■ M is called a local observable if there exists A G 
such that / is J^AxG-measurable, i.e. if /(a) = f{ri) whenever cta = i/a- 

ii. An event A G is called local if is a local observable, i.e. if A G 1Jagr° ^AxG- 

iii. An event A G is called a-local if it is the countable union of local events. 

2.1.3. Topological structure. Physically, two configurations are close whenever inside some 
large compact set each configuration is a slight deformation of the other. This means that 
each particle inside this compact set of one configuration can be matched to a neighboring 
particle of the other and vice versa. The precise definition is as follows. 

Definition 2.9. 

i. Given 6 > 0 and ^,ri E J\f{S x G) we say that ^ is 6-embedded in p if there exists 

an injective application p : [^] —)■ [rj\ such that d{TisxG{.lxii)iT^sxG{p{lxii))) < d 
for all £ ['C]) with ttsxg '■ (S' x G) x N ^ S x G the projection onto S x G. 

We denote it by ^ :<s rj. 

ii. Given a particle configuration ^ G J\f{S xG), a compact set JF C S' x G and 5 > 0, 
the (K, (j)-neighborhood of ^ is the set 

(Ok,s = {h e A/'(S' X G) : Si 5 V and rjx SiS 

The topology of the configuration space is the one defined by these neighborhoods. 

Definition 2.10. The vague topology on M(S x G) is the topology generated by the basis 

® ^ ^ Af{S X G),K C S X G compact and 5 > O}. 

A number of observations are in order. 

Remark 2.11. 

(a) J\f{S X G) admits a metric consistent with the vague topology, under which it is 
a separable metric space. It is also complete whenever G is compact. 

(b) The (T-algebra iF defined in fl2.2p is actually the Borel a-algebra corresponding to 
the vague topology on M(S x G). 

(c) The vague topology is usually defined as the one generated by the neighborhoods 

(0/i,...,/n,5 = {n e AA(^ X G) : -v(fi)\ <d , i = l,...,n} 

for 6 > 0 and fi, ■ ■ ■, fn continuous functions vanishing outside some compact set. 
This definition is equivalent to ours, but we will not use it in the sequel. 

2.1.4. Point processes on SxG. We call any random element of M(S x G) a point process 
on S' X G. Every point process X on S' x G is characterized by its distribution Px, which 
is a probability measure on J\f(S x G); the original measure space on which the process is 
defined plays no role. In the sequel, we will study convergence of point processes on S' x G 
and, in general, of probability measures on J\f(S x G). Besides the well-known notion of 
weak convergence, in our work we will also consider the notion of local convergence, which 
we define now. 


Definition 2.12. 
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i. A sequence (/in)nGN of probability measures on Af{S x G) converges locally to a 
probability measure jj, on J\f{S x G) if 

lim /i„(/) = /i(/) 

n—>-+cxD 

for every bounded local function / : A/'(S' x G) —)■ M. We denote this by /i fi. 

ii. Likewise, a sequence (A'„)ngN of point processes on S' x G converges locally to a 

point process A on S' x G if Px^ Px- We denote this by A„ X. 

In the general setting local observables need not be continuous, so that local and weak 
notions of convergence may not coincide in general. However, since uniformly continuous 
functions can be approximated arbitrarily well in the supremum norm by local observables, 
we have that local convergence is always stronger than weak convergence. Both notions 
coincide in fact whenever S and G are both countable and discrete. 


2.2. Gas models. Our models consist of three ingredients: (1) a configuration space, (2) 
an underlying “free” measure on the configuration space, and (3) a notion of interaction 
encoded in Hamiltonians. The first ingredient was the object of the preceding subsection; 
here we discuss the remaining two and give the general definition of gas model. 


2.2.1. Poisson processes on S x G and the free measure. In all gas models the underlying 
free measure will be a Poisson distribution on J\f{S x G) for some appropriate intensity. 
We give the definition of Poisson distribution below. 


Definition 2.13. Let z/ be a measure on [S x G,Bsxg) of 5'-locally finite allocation [i.e. 
such that z/(A x G) < +cxo for every A G 

i. The Poisson distribution with intensity v is the unique measure on M[S x G) 
which satisfies 


7r^({e e M{S X G) : i{Bf) = ki for all i 


i,...,n})= n 
2=1 




for all ki,... ,kn ^ No, disjoint Hi,..., e Bsxg n G N. 
ii. A point process X is called a Poisson process with intensity u if it is distributed 
according to tt^, i.e. for every finite collection of disjoint sets Hi,..., H„ G Bsxg 
the random variables A(Hi),... ,X{Bn) are independent and Poisson-distributed 
with respective means p(Hi), ..., z/(H„). 


2.2.2. Hamiltonians. In this article we aim at treating general gas models, including also 
contour ensembles. The latter are not amenable to a standard Gibbsian description. 
Rather, the energy cost of particle configurations in these models is described in terms of 
a Hamiltonian prescription, according to the following formal definition. 


Definition 2.14. 

i. A Hamiltonian prescription on J\f{S x G) is a family of measurable functions 

H = {Ha|^ : AA(A X G) ^ (- 00 , +oo] : A G , r; G x G)}. 

The function H^^ is called the local Hamiltonian on A with boundary condition r]. 

ii. Given a Hamiltonian prescription H, A & Bg and a configuration rj G x G), 

the energy leap in A relative to rj is the function : A x G —?■ (—C) 0 , -|-oo] 
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defined by 

{ HA\^{riAxG + K) - if < +oo 

+CX 0 otherwise . 


The local Hamiltonian HAjr) measnres the energy cost of inserting confignrations in A 
if snrronnded by the confignration r]. The energy leap AEA\rj represents the energy cost 
of placing an additional particle 7 ^ inside A when in the presence of the confignration rj. 
In particular, we may dehne the impact relation ^ on S' x G by the rule 

(2.3) % 3 A e and 7 e M{S x G) with AEA\rj{^x) 7 ^ AEA\rj+ 5 ^^ (jx)- 

If 7 y ^ jx we say that has an impact on 7 a,. This relation is not necessarily symmetric. 


2.2.3. Definition of gas model. Every gas model is defined by a pair (i/, El), where n is an 
intensity measure and H a Hamiltonian prescription. The former describes how particles 
would be distributed if there would be no interaction among them, and the latter specihes 
this interaction. Our aim of considering general point processes forces us to list a relatively 
long list of assumptions. As illustrated below, these assumptions are naturally satished 
by usual examples. 


Definition 2.15. A gas model on S' x G is a pair {n, H) verifying the following conditions; 

(1) S-locally finite allocation. For all A G the measure n satishes z/(A x G) < + 00 . 

( 2 ) Diluteness condition: HA\rj{0) < +00 for every A E Bg and 7 G Af{S x G). 

(3) Existence of infinite-volume energy leap function: The limit 

AEfi-fx) ■■= lim AEA|r,( 7 ^) 

A y' O 

exists for all G S' x G and 7 G J\f{S x G). 

(4) Bounded energy loss and allowance of particles: 


— oo< AE := inf 


agbo 


inf^ AEa\^{^x) 

Ax^AxG 

^■n^AfiSxG) 


< +CX0. 


(5) Integrable interaction range: Let the interaction range of B E BgxG be the set 


I{B) = { 7 y GS'xGidya-Gi? such that fiy . 

Then I{B) is measurable and i^(/(A x G)) < +cxd for each A G 

(6) Measurability of local Hamiltonians: Given A E Bg and y^, G S' x G, 

i. The application (^, 7 ) LfA|7?(0 is (JAxg G) J'(A'=xG)n 7 (AxG))-measurable. 

ii. The application rj i-A AEyipfx) is J3'(7,„)-measurable. 


These conditions are satisfied by all physical systems of interest we know of (although 
for some discrete systems, like the Ising model for example, one may need to consider an 
alternative lattice gas representation for these to hold), with the exception of condition (4). 
This condition is violated, for example, by interactions of Lennard-Jones type. Indeed, 
if one considers a ring of particles of a radius for which the L-J potential is negative, 
then the addition of a particle at the center of the ring would lead to an energy leap that 
becomes arbitrarily low with the (potentially unbounded) number of particles in the ring, 
thus yielding AE = — 00 . This suggests that, except for systems with purely nonnegative 
interactions, the validity of the leftmost inequality in (4) is tantamount to the existence 
of some sort of hard-core requirement preventing arbitrarily large amounts of particles 
inside bounded regions. 
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2.2.4. Gas kernels. Every gas model defines a family of probability measures on J\f{S x G), 
called gas kernels, which describe the local behavior of the system in bounded volumes. 
We introduce this family of gas kernels below. 

Definition 2.16. The gas kernel of (z/, H) on the volume A ^ Bg with boundary condition 
7 ] G Af{S X G) is the probability measure /iai?? on Af{S x G) given by 

(2.4) k'Alri ^A\r) ^ 

where we make the identihcation J\f{S x G) = A/'(A x G) x x G) and ojf denotes 

the probability measure on A/{A x G) dehned through the relation 

g ^ A I rj 

dUAlr, = -7^ - dirf 

■^A\r] 

with TTf denoting the Poisson distribution on A/{A x G) with intensity paxg cind 

ZA\rj= [ 

JAfiAxG) 

being the normalization constant. Notice that by assumptions (l)-(2) in Dehnition 12.151 
we have 

Za\v > = 0) = e-^G^^^)+HAir,m) > g 

so that bJA\ri is well dehned. 

The measures /xait? describe the local behavior of the system inside the volume A when 
the conhguration outside A is hxed as rj. The true objects of interest for us are, however, 
the possible local limits of these along suitable boundary conditions. 

Definition 2.17. Let {v, H) be a gas model. 

i. We say that a conhguration rj G A/{S xG) has hnite TT-interaction range whenever 
r 7 (/(A X G)) < +00 for every A G Bg, i.e. if rj has only hnitely many particles 
interacting with those in any given bounded volume. 

ii. A probability measure /i on A/{S x G) is called a (z/, H)-gas measure if there exist 
T] G M{S X G) with hnite iL-interaction range and (A„)„gN ^ B^g with An S 
such that 

(2.5) fiA„\r, ^ /i. 

By condition (5) in Dehnition 12.151 conhgurations with an inhnite iL-interaction range 
are not physically admissible for the system. Thus, allowing these as boundary conditions 
in fl2.5p may lead to pathological limits of no physical meaning. Indeed, it is not hard to 
device examples in which <50, the ^-measure on the empty conhguration, can be obtained 
as the local limit in fl2.5p along a boundary condition with an inhnite TT-interaction range. 
This is why in the dehnition of gas measure we disregard this type of boundary conditions. 

2.2.5. Gibbsian gas models. Gas measures possess a clear physical interpretation for the 
particular case of Gibssian models. These models are dehned by Hamiltonian prescriptions 
in which local interactions do not depend on the particular volume under consideration. 

Definition 2.18. We say that a model (z/, H) is Gibbsian if the Hamiltonian prescription 
H satishes the consistency property 

(2.6) HA\ri{o') = LfA|fTA-»?A= (“^a) a .^(A\A)|0a'??ac (‘^AXa) 

for every A C A G Bg, a G Af{A x G) and rj G Af{S x G). 









STABILITY OF GAS MEASURES UNDER PERTURBATIONS AND DISCRETIZATIONS 


9 


Notice that condition fl2.6p implies that, as we anticipated, the energy leap AE^r]{lx) 
does not depend on A. In particular, assumption (3) in Dehnition 12.151 immediately holds 
for Hamiltonians verifying fl2.6p . Furthermore, condition fl2.6p yields gas kernels satisfying 
a consistency relation of the form 

(2-7) ^^A\r| — J f^A\^ 

for every pair of volumes A C A G and rj G J\f{SxG). A family of gas kernels satisfying 
(12.7p is called a specification. The inhnite-volume measures relevant to Gibbsian models 
are usually introduced through the A —)■ S' version of (12.7p . 

Definition 2.19. Let {u, H) be a Gibbsian model. A probability measure /i on J\f{S x G) 
is called a Gibbs measure for (z/, H) if for every A G 

( 2 . 8 ) = j fTAindfriri). 

Thus, Gibbs measures are precisely those inhnite-volume measures which are consistent 
(in the sense of (I2.7p l with the local description of the model given by its specihcation. 
Hence, we may think of them as the measures describing the global states of our system. 
The relation between Dehnitions 12.191 and 12.171 follows from the fact that, for a large class 
of Gibssian models (see |21| for details), any local limit of gas kernels as put in (12.Sp is 
in fact a Gibbs measure for {u, H) and, furthermore, all extremal Gibbs measures can be 
obtained in this way. 

In practice, almost every model of physical interest is Gibbsian. There is, however, 
one important exception: contour ensembles. These are of particular relevance since they 
constitute one of the main tools for studying Gibbsian systems in the low-temperature (or 
condensation) regime. We give further discussion on Gibbsian and non-Gibbsian models 
in the following section. 

2.3. Examples of gas models. We now present some examples of models to illustrate 
the dehnitions of the previous section. Later in Section 0] we will also use these as ground 
for applications of our results. 


2.3.1. Models given by an interaction potential. The typical way in which Hamiltonians 
satisfying (12.6p are specihed is via an interaction potential, i.e. a family $ = of 

symmetric functions : (S x G)” {—oo, -|-cxo] subject to appropriate measurability 

and summability requirements so that the local Hamiltonians 


(2.9) 


H. 


Mv 


a := 


E 


n\m\ 


n>l 

Tn>0 


^ ^(n+m)/ (1) (n) ~(1) ~(m)\ 

/ ^ ^ \ Ix ^ ^ fx ^ fy ^ ^ fy ) 


are well-dehned and satisfy all the pertinent conditions in Dehnition 12.151 The resulting 
Hamiltonian prescription is said to be specified by <F and will be often denoted by H^. 
The function is called the n-body interaction of the potential <h. 

We say that a model (z/, H) is given by an interaction potential ii H = for some <F. 
It follows from (12. 9 p that any such model is Gibbsian. Furthermore, we have that; 

• ^ ^ ^ h ^ d\f{S X G). 

• For any G S' x G and rj G Af{S x G) the energy leap takes the form 


( 2 . 10 ) 


ae: 


'hJ = E — 

ml 


E 




a(i) 

’ D ’ 


• • • , 7 , 


(m) 


)■ 


m>Q 
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• The bounded energy loss condition in (4) of Definition 12.151 in this case reduces to 
the existence of a constant C > 0 such that 






^( 1 ) 
X) Jy ) 


’ D 


)>-C 


,) 6 [; 


for every G S'x G and rj G J\f{S x G). This condition is well-known and standard 
in the study of gas systems (see [201 Section 1.2]). 

Models specified by an interaction potential are the most common among gas models. 
Below we give some examples. 


2.3.2. The discrete Widom-Rowlinson model. It is a classical hardcore interaction model, 
first introduced by Lebowitz and Gallavotti in [16]. It involves particles of two types, say 
(-l-)-particles and (—)-particles, located at the sites of the discrete lattice 1/ for d > 1. 
The interaction between particles allows at most one particle site and forbids any two 
particles of different type from being within a certain fixed distance fc G N of each other. 
The corresponding gas model is defined by the following ingredients: 

• Location space S' = and spin set G = {+, —}. 

• Intensity measure given by 


( 2 . 11 ) 


u — A_|_ * Cjd X (5_|_ “t“ A_ * C'^d X , 


where A+, A_ > 0 are two fixed parameters known as the fugacities of (±)-particIes 
respectively and c^d denotes the counting measure on Z'^. 

• Hamiltonian prescription H specified by the potential given by 


( 2 . 12 ) 


<*®(7i,7») : = 


-|-cx) ii X = y 

-|-cx) if 0 < ||x — |/|loo < k and 7 7 ^ 7 
0 otherwise. 


Alternatively, one could define the model by considering instead the pair (c 2 dx{+,-}, Lf), 
where Cgdx{+^_} is the counting measure on Z'^ x —} and H is specified by the potential 
$ = $(2)) with as in fl2.12p and 


[ - log A+ if 7 = 

•i“>(7j = < 

— log A_ if 7 = —. 

Both representations are equivalent in the sense that they produce the same gas kernels. 
Nevertheless, for our analysis it will be more convenient to adopt the first representation. 
The reason for this choice will be explained later in Section HI We adopt this representation 
also in the remaining examples. 


2.3.3. The continuum Widom-Rowlinson model. In the continuum version of the model, 
particles are now located throughout the entire Euclidean space and the interaction 
forbids particles of different type from being within a certain distance r > 0 of each other. 
It was originally introduced by Widom and Rowlinson in [21| and later studied in mm- 
Its formal ingredients are: 

• Location space S' = and spin set G = 
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• Intensity measure 

1 / := A+ ■ X (5+ + A_ ■ X 5_, 
where is the Lebesgue measure on 

• Hamiltonian prescription H specihed by the potential $ = given by 

{ +CX 0 if Ijx — y\\oo < r and 7 7^ 7 

0 otherwise. 

Notice that the hrst term in fl2.12p excluding multiple particles in one site is now missing 
from fl2.13p . This is because the Poisson distribution already assigns zero probability to 
conhgurations with more than one particle per site and so this term becomes unnecessary. 


2.3.4. The Widom-Rawlins on model with generalized interactions. Several generalizations 
of the Widom-Rowlinson model are worth looking into. One interesting possibility is to 
consider a model in which nearby pairs of particles of opposite type are not necessarily 
forbidden, but merely discouraged, and also intra-species repulsion terms are included. 
Such generalization is dehned through decreasing functions h,j± : M’*' —)> [0, -|-cxd] with 
bounded support, by replacing (I2.12p - fl2.13p with the 2-body interaction 




h{\\x-y\\oc) if 7 ^7 

J-dIa:-i/lloo) if7 = 7 = - 
J+(ll^-l/lloo) if7 = 7 = +. 


We call (h, j_, j+) the repulsion vector. The original continuum Widom-Rowlinson model 
is obtained by setting h := (-|-cxD)l[o^r] and j± = 0. The discrete version corresponds to the 
same choice of h and j± := (-|-oo)l{o}. We refer to [9] where these type of generalizations 
were investigated. 


2.3.5. The thin rods model in 1?. Given fc G N we consider a system of hard rods in 
of zero width and length 2k whose centers are located at the sites of 1?. Each rod has an 
orientation specihed by an angle 7 G [ 0 , vr) with respect to the x-axis, and the interaction 
forbids any two rods to intersect. More precisely, if for r > 0 we set 


:= {f ■ (cos7,sin7) : f G [-r, r]} 


then the thin rods model in T? is dehned by: 

• Location space S = T? and spin set G = [0, tt). 

• Intensity measure := A • x p, where A > 0 is called the fugacity of rods and 
p is a probability measure on G called the orientation measure. 

• Hamiltonian prescription H specihed by the potential <I> = given by 


(2.14) 




-Fco if (L^ + x) n {L^ + y) ^ (I) 
0 otherwise. 


Of particular interest to us is the case when the orientation measure is given by 


(2.15) p = pSq + {1-p)S- 

for some p G (0,1). This model is identical to the discrete Widom-Rowlinson model, with 
the exception of an additional repulsion term between particles of the same type. Indeed, 
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by identifying G with {+.—} we have that u equals fl2.1ip for A_|_ := pX and A_ := (1 —p)X, 
while in (I2.14p can be rewritten as + where is as in (12.12p and 




-|-oo 

if 7 = 7 = -|- and \xi — yi\ 

< k 

(2.16) 

= < 

-1-00 

if 7 = 7 = — and \x 2 — 1/2 

< k 



0 

otherwise. 



We call this particular system the nematic thin rods model. More details about this model 
can be found in min] and references therein. 


2.3.6. The thin rods model in Similar to the previous model, the rod centers are now 
located at arbitrary points of and the rod lengths are 2r for some fixed r > 0 which is 
not necessarily an integer. The model is formally dehned by; 

• Location space S' = and spin set G = [0, tt). 

• Intensity measure u = X ■ x p, where A > 0 and p is the orientation measure. 

• Hamiltonian prescription H specihed by the potential $ = where 


:= 


+00 if (L!^ + x) n (Lt + I/) 0 

0 otherwise. 


For p as in fl2.15p . the model is equivalent to the continuum Widom-Rowlinson model since 
TT^ assigns zero probability to conhgurations for which in (I2.16p would be nonzero. 
Thus, both the Widom-Rowlinson model and the nematic thin rods model have the same 
continuum version. We refer to |2| where the model with a hnite number of orientations 
was studied. 


2.3.7. The Peierls contours model. To conclude, we show an example of a model which 
is not Gibbsian, but still hts into our framework; the Peierls contours model. It was hrst 
presented by Peierls in [18] to study the Ising model at low temperature (see also |12|). 
For simplicity, we focus only in the 2-dimensional case. 

The Ising model is the lattice system on the conhguration space { — 1, +1}^ dehned by 
the set of hnite-volume specihcations ; A G ^^ 2 ,^ G { — 1, +1}^^} given by 

^A\r, 

where 

(2-17) ;=-^ a{x)a{y)-P Y 

x,y£A x(iA,y^A 

lh-2/l|2=l \\x-y\\2 = l 

for a hxed parameter /3 > 0 known as the inverse temperature. Let us observe that, since 
the interaction has range one, the boundary condition p is involved in fl2.17p only through 
its values on the external boundary dA of A dehned as 


dA := {yiA^-. d2{y,A) = 1}. 


Of particular interest are the boundary conditions -|- and —, corresponding to p{x) = -|-1 
and p{x) = —1 for all x G respectively. It can be seen that the local limits 


/i+ ;= lim Pai 

A /' Z 2 ^'■1 


p ;= lim Pai 

A /' Z 2 ^'■1 


and 
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both exist and constitute the unique extremal Gibbs measures of the model for a fixed /9, 
in the sense of Definition 12.191 fsee |10) and references therein for details). If [3 is such that 
/i"*" and do not coincide, we say that a phase transition occurs at inverse temperature (3. 
Peierls showed the existence of a phase transition for all sufficiently large values of (3 by 
considering the following geometric description of configurations in terms of contours. 

We begin by fixing + as the boundary condition and letting A G B ^2 be a square. 
Now, consider := Z^ + (|,|), the dual lattice of Z 2 . Given an edge e joining two 
neighboring sites in Z^, let e* denote the unique edge joining neighboring sites in Z^ 
which is orthogonal to e. We call e* the dual edge of e. Furthermore, consider: 

• e(A), the set of edges in with at least one endpoint in A. 

• e*(A) := {e* : e G e(A)}, the set of dual edges of e(A). 

• A*, the set of sites in Z^ which are endpoints of edges in e*(A). 

Given a configuration a satisfying the boundary condition + outside A, let denote 
the set of dual edges e* G e*(A) such that e joins two sites x,y with different spin, i.e. 
a{x)a{y) = —1. With a little work it is possible to show that the edges in Do- join up to 
form closed curves (which may contain loops). This set of curves can be decomposed into 
connected components 71 ,..., 7 ^. We call any of these components 7 * a contour^ and write 
Fo- ;= { 71 ,..., 7 n} for the set of contours of a. It can be seen that the assignation a 1 —)■ Fg- 
is in fact a bijection: given a finite family F of mutually disjoint contours contained in 
e*(A), there exists a unique configuration crp satisfying the boundary condition + in A'^ 
which has F as its set of contours. Furthermore, if I 7 I denotes the number of edges in 7 , 
then for any such a we have 


( 2 . 18 ) 



_e-2/3E^6r<, bl 

IFa 


where Wa is a normalizing constant depending solely on A. Thus, whenever A is a square, 
with fl2.18p we obtain an alternative representation of P^a\+ terms of a system of contours 
interacting by exclusion. In the current framework of gas models, this system is defined 
by setting: 

• The dual lattice Z^ as the location space S. 

• The set of contours rooted at the origin 0* := (|, |) G Z^ as the spin set G. Here, 
we say that a contour 7 is rooted at a: G Z^ if a; is the smallest site belonging to 7 
(with respect to the lexicographical order). Thus, we interpret any 73 . G Z^ x G 
as the contour shape 7 rooted at x. 

• The intensity measure u given for each 7 ^; G S' x G by z/( 7 a;) := 

• For each A G and rj G A/'(Z^ x G), the Hamiltonian specified as in (12.9h 

but for the local potential <hA = (*^a\ $( 2 )^ given by 

( +CX) if 7 a; n 7 y 7 ^ 0 

$^"H7x,7x) = < 

I 0 otherwise 

and 

( +CX) if 73 ; n (Z^ - A) 7 ^ 0 

<1>aGJ = I 

I 0 otherwise. 

The interaction term is responsible for the exclusion among different contours, 
while the term bans those contours which are not contained in A. 









14 


ROBERTO FERNANDEZ, PABLO GROISMAN AND SANTIAGO SAGLIETTI 


The resulting pair (z/, H) is called the Peierls contours model. The main physical interest 
of this model lies in the fact that phase transitions in the Ising model at low temperatures 
can be understood in terms of the diluteness properties of gas measures for {v, H). Indeed, 
if : A G ^ is the family of kernels induced by (z/, H), then fl2.18p 

can be rewritten as 


(2-19) hA|+M = AL|0(r-) 

for any square A G B ^2 and spin configuration a G { —equal to + outside A. 
Here, 0 denotes the empty contour conhguration. Using fl2.19p and the spin-flip symmetry 
of the Ising model, one can show that if the gas measure jj,^ := liniA/iza is sufficiently 
diluted (which occurs at low temperatures) then the inhnite-volume measures /i^ and 
in the Ising model are distinct. Therefore, by changing the (local) spin variables into new 
(non-local) contour variables, the proof of the existence of a phase transition reduces to 
the proof of some form of diluteness of the contour measure. We refer to [21] for details. 

In our present context, however, this model is of interest also for another reason; it 
constitutes the canonical example of a (physically relevant) non-Gibbsian system. Indeed, 
for any G x G, A G B ^2 and y G A/'(Z^ x G) we have that the energy leap 

( -Fcx) if 7 a; n (y) 7 ^ 0 or fl (Z^ - A) 7 ^ 0 
AU'A|r)(ya;) = A 

I 0 otherwise 


depends on the volume A through the restriction imposed by the interaction term 
implying that the model is not Gibbsian. However, the limit 


lim AUA|^(ya;) 

V/'Z2 '' 


-Fcx) if yo; n (y) 7 ^ 0 
0 otherwise 


exists for all choices of rj and 73 ., so that the conditions in Dehnition l2.15l are still satished. 
The restriction imposed by ^ must be, nonetheless, included for (I2.19P to hold, as 
is supported on conhgurations a such that Tg. is always contained in A*. 


2.3.8. Other examples. More in general, our treatment is also adapted to handle systems 
in the following general classes: 


General contour ensembles. The Peierls contours discussed above are particularly simple 
because the Ising Hamiltonian is symmetric under the overall flipping of conhgurations. 
More general non-symmetric cases are the object of study of Pirogov-Sinai theory |23] . 
The main features of contours dehned in this theory are the following: (i) Gontours are 
“thick" subsets formed by collections of plaquettes, (ii) contours include some additional 
information (color, conhgurations on both sides or, in general, the conhguration on the 
relevant plaquettes), (hi) each reference conhguration has a specihc contour ensemble, (iv) 
contour ensembles are not of physical nature and only external contours coincide with the 
physical “defects” in the presence of reference conhgurations, (v) contour weights include 
ratios of partition functions that must be bounded so to obtain exponential expressions 
similar to the one in fl2.18p . The last property is encoded in the expression “contours must 
satisfy a Peierls condition”. If these conditions are met, by proceeding as for the Peierls 
contours one can show that the diluteness of a contour ensemble implies the existence of 
a measure “tilted" towards the corresponding boundary conhguration. We observe that 
PS contours ensembles do not ht the Gibbsian framework for the same reason that in 












STABILITY OF GAS MEASURES UNDER PERTURBATIONS AND DISCRETIZATIONS 


15 


the Peierls contours model. Our results applied to these PS contour ensembles can yield 
not only proofs of the existence of phase transitions in the associated spin systems, but 
also the stability of the resulting phases with respect to perturbations and discretizations. 
This will be exploited in a subsequent publication j^. 

General polymer models. These models —introduced by Gruber and Kunz [13]— involve 
general geometrical objects subject to a general hard-core condition dehned in terms of a 
“compatibility" relation (see dais])- They are the traditional target of cluster-expansion 
or closely related methods [1]. Our results extend uniqueness and mixing properties to 
a larger region of parameters than expansion-based treatments, at the cost of sacrificing 
analyticity considerations. 

General point processes. Point processes are the genesis of the ancestor algorithm exploited 
in this paper. All the models presented in one of the original publications [ 8 | fall within 
the scope of our treatment: area-interacting processes, Strauss process, loss networks, 
random cluster model. In fact, our results apply to models that combine the generality 
of polymer models —that do not require a geometric underlying space— with that of 
point processes —that allow soft as well as hard-core interactions. A full presentation of 
the ancestor algorithm in such a general framework, and involving even weaker dilution 
requirements, is the object of a separate paper | 22 ) . 

2.4. Approximation families. We now describe the perturbations of the conhguration 
space that will be considered on the target model. These include, but are not limited to, 
discretization schemes on both the location space and spin set. 

Definition 2.20. A family D = {Ds)e>o of measurable applications : S x G ^ S x G 
is called an approximation family if the following conditions are satished: 

i. For any B G and 5 > 0 there exists G Uo <£<5 ^ ■ 

11 . For any B G ^^sy.G ^ ^ H there exists B^^'j G ^^sy.G with Uo<e<5^e(5) C %)• 
iii. There exists a : M>o —)■ IF>o with hm£_,.o+ a{e) = a(0) = 0 such that 

c^ 5 xG(^£( 7 x), 7 a:) < a{e) 
for every E S x G and £ > 0. 

The application is called the e-approximation operator. 

We note that conditions (i)-(ii) are merely technical requirements needed for the proofs, 
the essence of Dehnition l2.2Ul is contained in (iii). In fact, if the metric structure on S' x G 
is such that Bs = { 73 , E S x G : d{jx, B) < 6} E Bg^Q holds for any B E Bg^Q then 
conditions (i)-(ii) are immediately satished and they can be removed from Dehnition 12.201 
On a side note, we observe that by dehnition Do is always the identity operator on S' x G. 
This is so for notational convenience. 

Examples 2.21. Some natural examples of approximation families include: 

• Spatial translations: Dehned on S' = given for each e > 0 by 

Deix,-f) = {x + e-v,-f) 
for some hxed unit vector v E M'^. 

• Spatial discretizations: Dehned on S' = M'^, given for each e: > 0 by 

= (x£,7) 
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where, for x = (xi,..., Xd) G we write 


( 2 . 20 ) 



• Spin rotations: Defined on G = the nnit sphere in d-diniensions, given for 

each e > 0 

DTix,-f) = (x,£-i?7) 
where R is some hxed rotation. 

• Spin discretizations: Defined on G = [0, tt), given for each e > 0 by 


( 2 . 21 ) 


Df(x, 7 ) = (x,e 



Approximations can be composed giving rise, for instance, to operators of the form 
that discretize both space and spin and may depend on more than one parameter. 
To simplify the exposition we will always assume that £ G M>o, but we point out that the 
extension to the case in which £ = (^i,..., Sk) is a vector of parameters is straightforward. 

In the sequel, for notational convenience we shall write 7 ^ instead of Moreover, 

for ^ G AA(S' X G) and each e > 0 we dehne G ATi^S x G) as the particle conhguration 

given by the standard representation 


^^(0 = X] 

7xG{«) 


The fact that is indeed locally hnite follows from (i) in Dehnition l2.20[ Furthermore, 

it follows from Lemma f6.3l below that hm£^o+ ^ vaguely. To simplify notation, 

we may sometimes write fy instead of 


3. Main results 


The main results featured in this article concern the particular class of heavily diluted 
gas models, which we introduce now. 


Definition 3.1. A model {n,H) satisfying Dehnition 12.151 is said to be heavily diluted if 
there exists a measurable function g : S' x G —)■ [1, +cxo) such that 

,-AE r 


(3,1) 


:= sup 


Ai.lx) 


R{ix) 




< 1 . 


The quantity 0 ^’^ is called the q-diluteness coefficient and q is called the size function. 
Heavily diluted models have a unique gas measure, as the following result in [2T] shows. 


Theorem 3.2 (Uniqueness of the gas measure in heavily diluted models). 

Let {u, H) be a heavily diluted model on S' x G. Then; 

i. The local limit /i := limA ^5 /iA|,, exists and coincides for any rj G J\f{S x G) with 
hnite iL-interaction range. In particular, (z/, H) admits a unique gas measure. 

ii. If (z/, H) is Gibbsian then p constitutes the unique Gibbs measure of the model. 


To state our results we need to introduce hrst the notion of negligible event on J\f{SxG). 

Definition 3.3. We say that an event N C J\f[S x G) is dynamically negligible for a given 
intensity measure z/ on S' x G if it satishes the following properties: 

i. is a (T-local Tr^'-null event. 
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ii. N is closed by addition of particles, i.e. rj ^ N,ri ^ ^ ^ ^ N, where rj ^ ^ 

whenever their standard representations satisfy {rj) C (^) and mri{'yx) < 
for every ^ Qrj- 

Examples of dynamically negligible sets will be given in the applications of Section 01 
As illustrated in Section 01 most realistic continuum models will not satisfy the limit (13.21) 
in the hypotheses of Theorem 13.41 below for every conhguration, but will rather do so only 
for conhgurations outside a dynamically negligible set. Thus, for our results to be of any 
real use in the continuum, it is necessary to allow for violations of (13.2p inside such sets. 
Fortunately, our results will still hold under this weaker hypothesis. Later in Section [6l 
we will introduce the even weaker notion of dynamically impossible sets when discussing 
further relaxations to the hypotheses of Theorem 13.41 below. 

We are now ready to present our main result. In its statement we add superscripts to 
the usual notation in order to identify the model which we are referring to. 


Theorem 3.4. Let (T*e)e>o be an approximation family and suppose that is 

a family of gas models such that: 

i. For every e > 0 the intensity measure is given by 

._ ,,o ^ n-i 


■= o D-\ 


ii. There exists a dynamically negligible set N for such that 


(3,2) 


lim AM-(73 = AM" (73 


ifO, 


e^ 0 + 


for every G S' x G and all p G J\f{S x G) with rj + 6^^ G W 


111 . 


AE := inf 

£>0 


■ 


■ 

inf 

agM 

inf AE^^eiAx) 

'yxEAxG ' ' 


o 

_7?gM(SxG) 



> —CX). 


IV. 


(3.3) 

(3.4) 

Then: 

(a) 

(b) 


There exists for each G S' x G a set Vi^jx) ^ ^SxG such that: 

• For every e > 0 one has the inclusion 

d ;' (/"'(73) c r(73. 

• There exists a size function g : S' x G —>■ [1, Tex) which verifies 


a 


u°,V 

:= sup 

7xe5xG 


-,-AE 


.Qhx 


'U(7x) 


g(yy)dz/°(% 


< 1 . 


Each model admits exactly one gas measure /i^. 

As e: —)■ O’*", we have the weak convergence 


e 0 

/i -)■ /i . 


(3.5) 


(c) There exists a coupling of the measures such that for any B E B[ 

there exists (a random) > 0 verifying that for all e < Eb 

Z% = DA zA 


0 

SxG 


D7Hb)J ■ 


In particular, Z^ Z^ with respect to the vague topology. 
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We point out that, in all common situations, the condition < 1 alone is enough 

to guarantee the validity of (iv) in the statement of Theorem 13.41 Indeed, as we shall see 
in Section 01 a set V(j^) satisfying fl3.3p can generally be obtained by slightly enlarging 
{'jx) ill some appropriate manner. If is a continuous function of the parameters 

of the model and all the remaining models (z/^, are “sufficiently close” to H^), 

then performing this slight enlargement will yield a coefficient very close to 0 ^°’^°, 

so that fl3.4p holds. 

In principle. Theorem 13.41 deals only with perturbations of the intensity measure which 
are given by approximations in the sense of Dehnition l2.20[ However, one may cover other 
cases of interest as well by hrst transferring perturbations in the intensity measure to an 
effective Hamiltonian prescription and then applying Theorem 13.41 In this way, we obtain 
the following important corollary, dealing with absolutely continuous modihcations to the 
intensity measure. This scenario typically represents perturbations in the parameters of 
the model: fugacity of particles, inverse temperature and interaction range among others. 


Corollary 3.5. Let (z/^,Tr^)e>o be a family of diluted models such that: 

i. There exists an intensity measure u on S x G such that v for every e > 0. 

ii. There exists a dynamically negligible set N for v such that 

lim AE^( 7 ^) = AEJ( 7 ^), 

£->■ 0 + 

for every 'jx ^ S x G and all rj G N{S x G) with r] + 5^^ G iV^, where 


AE^( 7 „) := |im AEX|^( 7 x) 

with 

^^A|r,(7x) := AE^^( 7 ^) - log ■ 

hi. 





■ 

AE := inf 

inf 

inf AE^(73,) 

'yx^^^G 


£>0 




0 

_r)GAf{SxG) 



iv. There exists for each 73. G S' x G a set Vi^jx) ^ ^SxG such that: 
• For every £ > 0 one has the inclusion 

(3.6) I^\<)^V{^x). 


• There exists a size function g : S' x G —>■ [1, + 00 ) which verifies 


vy 

a ’ := sup 






< 1 . 


Then: 

(a) Each model admits exactly one gas measure /i^. 

(b) As e —)■ O’*', we have the local convergence 




e 





(c) There exists a coupling of the measures such that for any B G 

there exists (a random) > 0 verifying = Z^ for all e < £ 5 . 
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4. Applications 


We discuss here some consequences of our main results for the models in Section 12.31 
First, we will focus on applications of Corollary 13.51 and then consider other applications 
using Theorem 13.41 in its full generality. The reader should keep in mind that, although 
we do not state it explicitly in each application, for every local and/or weak convergence 
of probability measures throughout this section Theorem 13.41 guarantees the existence of 
a coupling in which the convergence takes place almost surely. 

4.1. Applications of Corollary 13.51 We illustrate the use of Corollary 13.51 bv showing 
the continuity in the parameters of the unique gas measure for the Widom-Rowlinson and 
Peierls contour models dehned in Section [T3l 


4.1.1. The Widom-Rowlison model. 


Theorem 4.1 (The discrete Widom-Rowlinson model). Given k eN and > 0, 

consider the discrete Widom-Rowlinson model with fugacities A^ and exclusion radius k. 
Then, if 

(4.1) «wkA(/, Aq , k) ;= max{Ao , A(/}(2/c + 1)“^ + min{Ao , A^} < 1 

there exists an open neighborhood U of (A"*", A“) such that for any (A’*', A“) G U the model 
with fugacities A^ and exclusion radius k admits a unique Gibbs measure 
Moreover, the application (A"*", A“) eA X~) is continuous on U in the local topology, 

i.e. for any (A/", A“) G U the following local limit holds: 

lim _ /i(A+, A") =/i(A+, A/). 

(A+ ,A )—^(A,,, ,A* ) 

Proof. Observe that for any 7 ^^ G x {+, —} and 7 G x {+,—}) we have 


(4,2) 

lirtx) = {lx} u {% : 7 = 

V 

VI 

1 

u- 

1 

and 




r -Loo 

if Tl{I{lx)) > 0 

(4,3) 

^Enilx) = { 



1 0 

otherwise. 


It follows from this that for the model , H) with fugacities A^ and exclusion radius k 

the associated diluteness coefficient in fl3.ip for any constant size function q is 

'^ = sup ( 1 ( 73 ;)) = max{A“, A’''}(2fc -I- 1)'^ -I- min{A“, A’*'}. 

'yx€Z‘^x{+-} 

Therefore, if a^j^{XQ , Xq , k) < 1 then there exists an open neighborhood U of (A’'',A“) 
such that ’^ < 1 for all (A''’,A“) G U. In particular, by Theorem 13.21 there exists 

a unique Gibbs measure X~) of the model (z/^^’^ , H) for any (A^, A“) G U. 

To see that (A"*", A“) i-A ^{X'^, X~) is locally continuous in U, we £x (A+, A^) G U and 
check that for any sequence (A+, A“)„gH ^ U converging to (A+, A^) we have 


MA+ A-)^MAi„A-). 

For this it suffices to see that , ih)„gNu{oo} satisfies the hypotheses of Corollary 13.51 

But notice that if for n G NU{oo} we write n"' := then -C with density 

given by 


dz/" 


(7x) 




}• 
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In particular, we have that for every G x {+, —} and rj G x {+, —}) 


(jjA' 

AE;(7j=A<(7j-log(^(7. 


du'^ 


A<(7.) = Ai?-(7j 


and also that 


AE = - log ( sup [max {A+ - A+, A„ - A^}] 

VnGN 


> —CX3 


if (K, K) is sufficiently close to (A+, X^). Furthermore, if for each G x {+, —} 
we choose V (y^) := /^(y^) then the inclusion fl3.6p immediately holds for all n G NU{cxo} 
and 




provided that (A+, A„ )neN is sufficiently close to (A+, A^^) so that AE is close enough to 0. 
Hence, we see that the sequence ^ H)n^mj{oo} verifies the hypotheses of Corollary 13.51 

(with the limit n —)■ cxo replacing the usual e —)■ 0) and this concludes the proof. □ 


Remark 4.2. If one uses the alternative representation of the Widom-Rowlinson model 
given by the pair (i>, H) in Section Ti. 3. 21 then for any constant size function q one obtains 
the larger diluteness coefficient 

= max{A“, X'^}{{2k + 1)'^ + 1) 

which leads to a smaller uniqueness condition than the one in Theorem 14.11 This is why 
we chose to include the fugacities in the intensity instead of the Hamiltonian prescription. 


To illustrate the need to consider dynamically impossible events, let us treat the case of 
the continuum Widom-Rowlinson model. We have the following analogue of Theorem 14.11 

Theorem 4.3 (The continuum Widom-Rowlinson model). Given ro > 0 and A^j", Aq > 0, 
consider the continuum Widom-Rowlinson model of fugacities A^ and exclusion radius lq. 
Then, if 

(4.4) a[[)^^(A([, Ao,ro) := max{Ao , A(^}(2ro)'^ < 1 

there exists an open neighborhood U of (AQ',AQ,ro) such that for every (A’'',A“,r) G U 
the model of fugacities A^ and exclusion radius r has a unique Gibbs measure /i(A+, A“, r). 
Moreover, the application (A+,A“,r) i-A /i(A+,A“,r) is locally continuous on U. 

Proof. The first assertion follows as in the proof of Theorem 14.11 by noticing that now 
fl4.2p is replaced with 

HXx) = {% : 7 = -7 , Ik - y\\oo < r}. 

To see the local continuity, observe that if (A+, A“, rn)nm converges to (A+, A^, roo) G M>o 
and for each n G N U {cxd} we denote by H'^) the model with parameters (A^, A“, r„) 
then the convergence 

AE;( 7 .) = AEf ( 7 .) - log (^^(7x)j ^ A<°“( 7 .) = AE-{^^) 
may not hold if y -|- 5^^ is inside the set 

A = {^ G A/'(M'^ X -}) : 3 % 7 ^ 7 ^ G (0 such that y = -y and \\y - z\\oo = lq}. 

Indeed, the problem lies when y^, is at a distance ro from a particle in rj of opposite type. 
If this is the case and r„ ^ Voo then we can have AE^" {'jx) = 0 ^ -|-oo = (y^,). 

However, since A is a dynamically negligible event for z/°°, one may disregard these cases 
and still proceed as in the proof of Theorem 14.II to conclude the result. □ 
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Similarly, we can obtain an analogons resnlt for the model with generalized interactions. 
For simplicity, we only state it for its continnnm version and omit the details of the proof. 

Theorem 4.4 (The continnnm Widom-Rowlinson model with generalized interactions). 
For each e > 0 let ns consider the continnnm Widom-Rowlinson model with fngacities 
and repnlsion vector (h^, ji, j+). Assnme that the following conditions hold; 

i. lime^o+ A+ = A([ and hm£^o+ K = K- 

ii. hme_j,o+ h^{r) = h^{r) and lime_j,o+J±(r’) = j^(r) for £-almost every r > 0. 

hi. hme_^o+ and lim£^o+ where 

nihe = snp{r > 0 : h^{r) ^ 0} and rrij^ = snp{r > 0 : i%{r) ^ 0}. 

iv. h^,j-,j+) < 1, where for any A"*", A“ > 0 and repnlsion vector (h, j_, j+) 

we define 


q;vfr(A"'^, h, j) := 2'^max{A m'^_ + + X max{m^, 

Then for e > 0 snfficiently small there exists a nniqne Gibbs measnre of the associated 
continnnm Widom-Rowlinson model. Fnrthermore, we have the convergence ^ 

The nniqneness regions of parameters prescribed by the above bonnds are comparable 
—althongh a little weaker— to those obtained via disagreement percolation argnments 
(see, for instance, |10|h It is not clear, however, whether the latter method can lead to 
the stability resnlts bronght by onr approach. 


4.1.2. The Peierls contour model. 


Theorem 4.5 (The Peierls contonrs model). Dehne the coefficient 

f3p = inf{/3 > 0 : ap{l3) < 1}, 

where 


(4.5) 


ap{(3) := snp 

7:.6Z2xG 





Then, ii P > f3p the model at inverse temperatnre P admits a nniqne gas measnre /i(/3). 
Moreover, the application P i-A /i(/9) is locally continnons on (0,/9p). 


Proof. Notice that the coefficient ap{P) is exactly the dilnteness coefficient for the model 
at inverse temperatnre P associated to the size fnnction q{'^x) '■= \lx\- Hence, nniqneness 
for P > Pp follows from Theorem 13.21 For the continnity, we proceed as in the proof of 
Theorem 14.11 Thns, we fix P^a > Pp and a seqnence (/3n)neN ^ '^>i 3 p converging to Poo 
and show that the family Tf)nGNu{oo} satishes all the hypotheses of Corollary 13.51 
This can be done as in the proof of Theorem 14.11 bnt one has to be carefnl that 


inf 

neN 


■ 


^ A 

■ 

■ 

inf 

O'x ^SxG 

-log 1 


= inf 

nGN 

inf 2(/3„-/3 oo)|7x| 

'yx^S X G 


= —oo 


if Pn < Poo for some n G N, so that if one takes v in the statement of Corollary 13.51 as 
then the same argnment that in the proof of Theorem 14.11 does not go throngh this time. 
To solve this, we take u := with := inf„gj^u{oo} Pn and then proceed as before. □ 
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We point out that the analogous result also holds for the model in higher dimensions. 
The condition in Theorem|4.5|for the d-dimensional model (d > 2) is satisfied for all /3>0 
such that 

(4.6) 

l>2d 

where W denotes the number of contours of perimeter For the two-dimensional model, 
this condition coincides with the Peierls condition presented in usual textbooks treatments 
of the Peierls argument. In higher dimensions, however, the Peierls condition is weaker, 
basically corresponding to changing the factor I in the left-hand side by Also, 

additional arguments based on the Borel-Cantelli lemma show that it suffices to have 
“-|-oo” instead of “1” in the right-hand side of fl4.6p in order to guarantee a phase transition. 
Of course, our theorem implies further properties that cannot be directly obtained from 
the original Peierls argument. 

4.2. Applications of Theorem 13.41 We now discuss some further applications which 
use Theorem 13.41 in its full generality. We shall focus in discretization schemes applied to 
the thin-rods and Widom-Rowlinson models and determine diluteness regimes in which 
the Gibbs measure of the continuum model is the scaling limit of the equilibrium measure 
in the discretized model. 


4.2.1. The Widom-Rowlinson model and its “universality" class. 

Theorem 4.6. For Ao,ro > 0 such that Ao(2ro)'^ < 1 we have the following: 

i. The continuum Widom-Rowlinson model with fugacity Aq and exclusion radius tq 
in admits exactly one Gibbs measure, which we shall denote by /i°. 

ii. The discrete Widom-Rowlinson model with fugacity e'^Ao and exclusion radius ^ 
in admits exactly one Gibbs measure if 0 < e < (2'^ro)“s — To. 

iii. Provided that 0 < e < (2'^ro)“3 — tq, as e —)■ 0 we have 

~e —1 0 

W O 1 /i 

where for each e > 0 we define the shrinking map : Z'^ x {+,— } {+,—} 

by the formula 

G(a;,7) := 


Proof. Assertions (i)-(ii) follow at once from Theorem 13.21 and (I4.4p - (l4.ip . respectively. 
To show (iii), we consider the spatial discretization family (Zlg®)e>o given by fl2.20p and 
for each e > 0 set: 

• The intensity measure W as W := z/° o where u := AqT'^ x (5+ -|- (5_). 

• The Hamiltonian prescription as the one given by the 2-body interaction 


(4.7) 


:= < 


-|-oo 

-l-cxo 


0 


ii X = y and e > 0 

if 0 < llx - y\\oo < To and 7 7^ 7 

otherwise. 


Notice that is precisely the continuum Widom-Rowlinson model of fugacity Aq 

and exclusion radius vq, whereas for each e > 0 the pair constitutes essentially 
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the ^-shrunken discrete model of fugacity s’^Xq and exclusion radius —. More precisely, 
for every A G and e > 0 we have 

(4-8) ht(A)|0 = hA|0 ° ^ 

where /i^|0 is the Boltzmann-Gibbs distribution with empty boundary condition associated 
to the discrete Widom-Rowlinson model and is the one associated to H'^). Thus, 
by taking the limit as A ^ Theorem 13.21 yields for 0 < £ < (2'^ro) d - ro 

= ffo 

where is the unique Gibbs measure of the gas model given by the pair (z/^, H^). Hence, 
to conclude the desired convergence it suffices to show that the family is under 

the hypotheses of Theorem 13.41 But notice that: 

• (i) in the hypotheses holds trivially by the choice of measures 

• hm£^o+ for all rj G J\f(R‘^ x {+, -}) and G x {+, -} 

such that rj + 6^^ is outside the dynamically negligible set A^i U N 2 for where 

(4.9) A^i = {^ G Af{W^ X {+, —}) : (j{{y} x {+, —}) > 1 for some y G M'^} 

and 

N 2 = [ie X {+, -}) : 3 %, 7 ^ G (0 with 7 = -7 and \\y - z||oo = lq}. 

• /S.E := mi^>Q = 0 by the repulsive nature of the interaction. 

• If for 73. G X {+, —} we take 

ViXx) ■= {(2/, 7) : \\y-x\\oc < 5 }u{(z,-7) : ||z - x||oo < Lq + < 5 } 

with (5 > 0 sufficiently small so that for any constant size function q we have 
af= sup zx°(H( 7.)) = Ao((2(ro + 6 )Y + (2<5)") < 1, 

7r^6K‘^x{+,-} 

then YiY) ^ for every e > 0 sufficiently small. 

Thus, the hypotheses of Theorem 13.41 are satished and from this the result follows. □ 

Theorem 14.61 in particular shows that, for the case of the Widom-Rowlinson model, in 
order to simulate the Gibbs measure of the continuum system the method of sampling from 
the Gibbs measure of the discrete model obtained by discretizing both the interactions and 
the conhguration space indeed yields a faithful approximation of the desired distribution, 
at least whenever in the heavily diluted regime and the discretization is sufficiently rehned. 
It is clear that Theorem 13.41 guarantees that the same statement holds in general for other 
types of gas models as well. We point out that, in general, this procedure is not equivalent 
to directly discretizing the continuum measure. Indeed, the proof of Theorem 13.41 shows 
that fIb.Sp is in general as much as one can expect, although this is not enough to imply 
that for any hxed e > 0 the Gibbs measure of the discrete system coincides with 
the discretization of the continuum Gibbs measure. As a matter of fact, the latter may 
sometimes fail to be a Gibbsian at all, i.e. there exists no Gibbsian gas model for which 
is a Gibbs measure (see e.g. [16]). 

Another interesting interpretation the can be made of Theorem 14.61 is that it portrays 
the continuum Widom-Rowlinson model as a “universality class” to which many families 
of discrete models converge when properly rescaled. Indeed, the only relevant information 
about the discrete models used in the proof was that the additional interaction forbidding 
multiple particles per site vanishes as £ —)■ 0, which suggests that the same convergence 
should also hold for other discrete models with extra interactions which become negligible 
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in the continuum limit. At least for heavily diluted gas models, this is indeed the case and 
is essentially a consequence of Theorem 13.41 We illustrate this fact below by considering 
the particular case of the nematic thin rods model as an example. 


Theorem 4.7. For Ao,ro > 0 such that dAoUg < 1 we have the following: 

i. The continuum Widom-Rowlinson model with fugacity Aq and exclusion radius tq 
in admits exactly one Gibbs measure, which we shall denote by 

ii. The nematic thin rods model in with fugacity e'^Xq and rod length 2^ admits 

exactly one Gibbs measure if 0 < e < 2rQ + ^ — 2ro. 

iii. Provided that 0 <e<J2rl + ^^-2ro, as £ —0 we have 


/i" o y 1 




Proof. The proof is almost identical to that of Theorem 14.61 One only has to replace fl4.7p 
with 

^ +CX0 ife>0, 7 = 7 = + and |a;i — yi\ < tq 




+CX0 ife> 0 , 7 = 7 = — and \x2 — 1/2I < Lq 
+CX0 if 0 < ||x — y\\oo < tq and 7 7^ 7 


, 0 otherwise. 


and fl4.9p with N+ U A^_, where 

N+ = e X {+, -}) : 3 7j^ ^ 7^ G (0 with 7 = 7 = + and y^ = Zi} 

and 

N_ = X {+, -}) : 3 % 7^ 7^ e (0 with 7 = 7 = - and 1/2 = ^2}- 

The rest of the proof goes through exactly as before. We leave the details to the reader. 


□ 


4.2.2. The thin-rods model. For simplicity, we focus only on the continuum model in M^. 
We begin by establishing the continuity in the parameters of the model with a fixed hnite 
number of orientations. In the sequel, given a number /c G N, a vector of k orientations 
9 = {9i,... ,9k) G [0, tt)^ and a probability vector p = (pi,... ,pk) G we shall write 
pS,p for the orientation measure given by 

k 

i=l 


Theorem 4.8 (The thin rods model in M^). Fix /c G N and, given Ao,/o > 0, 6»o G [0,71)'= 
and a probability vector pq consider the thin rods model in with fugacity Aq > 0, 

rod length 21q > 0 and orientation measure . If 


(4.10) 


Iqi P 


,+ iPO 


) ;= 4A°(/°)" sup 

7€[0,7r) 


sin(7 — 7)|dp^°’^°(7) 


< 1 


then there exists an open neighborhood U of (Aq, /q) (^o,Po) such that for all (A, l,9,p) E U 
the model with parameters (A, I, 9,p) has a unique Gibbs measure /i(A, I, 9,p). Moreover, 
the applications {X,l,9,p) 1-7 p{X,l,9,p) and {X,l,p) 1—)■ p{X,l,9o,p) are respectively 
weakly continuous on U and locally continuous on = {{X,l,p) : {X,l,9o,p) E U}. 
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Proof. A straightforward calculation shows that if H^) is the model with parameters 
(Aq, Iq, 0o,Po) then the diluteness coefficient associated to any constant size function q is 


a, 


fi jjO _ 


= 4A°(/°)2 sup 


sin (7 — 


7€[0,7r) L'^O 

From this the first assertion about uniqueness of the Gibbs measure immediately follows. 
To establish the weak continuity, we fix {XooJoo,(^oo,Poo) ^ U and an arbitrary sequence 
(An,/n, 6 'n,Pn)nEN ^ U Converging to (Aoo,/oo, 6 * 00 ,^ 00 ) ^ u and show that 

/^(A/i,/ri; Pn) ^ /^(Acxd? ^oo: Poo\ 

For this purpose, we dehne the approximation family TX = (iA„)„eNu{oo} by the formula 

{x, On{i)) if 7 = 0oo{i) for some i = 1,..., k 


Dn{x,'j) = 




otherwise. 


where for each n G M U { 00 } and i = 1,... ,k, 6n{i) denotes the i-th coordinate of 6n{i)- 
Furthermore, for each n G N U {c)o} consider the Hamiltonian prescription given by 
the potential = (*hn \ 4)^ ^), where 


Ki'U.) = 



if 7 = 9n{i) for some i = 1,..., k 
otherwise. 


$(2)( ^if (g" + ^) (4" + y)^ 

^ ^ 1 0 otherwise. 


and 


Now, if we write := ^ where := Aoo>C^ x then it is not hard to see that 

the model is equivalent to the thin rods model with parameters (A„, G, 

The result now follows at once upon noticing that the sequence Ff"')neNu{oo} satishes 

the hypotheses of Theorem 13.41 if (A^, In, 9n,Pn)nm is sufficiently close to (Aoo, loo, 9oo,Poo)- 
Finally, to obtain the local continuity whenever the vector of orientations is fixed we follow 
a similar argument using Corollary 13.51 instead. □ 


To conclude, we consider the case of a continuum of possible orientations and derive 
the analogue of Theorem 14.61 in this context for spin discretizations. For dehniteness, 
we state the result for the continuum model with an uniform orientation measure on [ 0 , tt). 

Theorem 4.9. Given A, I > 0 , let us consider the thin rods model on with fugacity A, 
rod length 21 and uniform orientation measure pu on [ 0 , 7 r). Then, if 

aPj^{\,l,Pu) = -A/^ < 1 

TT 

there exists a unique Gibbs measure p of the model and for every n G N sufficiently large 
the model with orientation measure 


^ n—l 


Pu ■ 

n 


i=0 


also has a unique Gibbs measure /i” which, furthermore, satishes /r"' -4 p. 

We omit the proof of this result since it goes very much along the lines of Theorem 14.61 
but using the spin discretization family introduced in fl 2 . 2 ip instead. 
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5. The Fernandez-Ferrari-Garcia dynamics 


In this section we study the Fernandez-Ferrari-Garcia dynamics first introduced in [6]. 
In their work the authors focus on the Peierls contours model of Section 12.31 and show 
that, for a sufficiently large value of the inverse temperature /3, the unique gas measure 
of this contour model can be realized as the unique invariant measure of this dynamics. 
In [2T] it is shown that this result can be extended to the broader class of gas models, 
where the sufficiently large inverse temperature condition is replaced by the one in fl3.ip . 
This extension of the dynamics will be the main tool used in the proof of Theorem 13.41 
We give a brief overview of it now. 

Given a gas model (z/, H), we consider the following dynamics on N'{S x G)\ 

• At rate the birth of new animals is proposed with intensity given by z/. 

• Each proposed for birth will be effectively born with probability 

where rj is the state of the system at the time in which the birth of is proposed. 

• Every animal which has effectively been born will have an independent lifetime, 
with exponential distribution of mean 1. 

• After its lifetime has expired, each animal dies and vanishes from the configuration. 
Thus, the infinitesimal generator L associated to such dynamics is formally given by 

(5-1) L(/)(a) = j + f{-ix))dp(ax) 

lx€{cr) 

for any bounded local function / : J\f{S x G) —?■ M. The following result is shown in |21| . 

Theorem 5.1. If [u, H) is heavily diluted then the Markov process with generator L exists 
and has the unique gas measure of (z^, H) as its unique invariant measure. Furthermore, 
there exists a time-stationary construction K, = of the process. 

For our purposes it will be convenient to understand the basics of the proof of this result, 
so we summarize it below. The idea is to construct /C explicitly as a suitable thinning of 
a non-interacting birth-and-death process on S' x G with the appropriate rates. 

To make matters more precise, let us consider the product space C = (S' x G) x M x M’*'. 
The elements of C are called cylinders, since any (73,, t,l) E C can be seen as a cylinder 
on S' X M of axis {x} x [t,t + 1] and diameter 7. However, we shall prefer to describe each 
cylinder G = s) G C in terms of its basis jx, its time of birth t and its lifespan 1. 

We denote these three features of G by basisiC) ,bc and Ic, respectively. 

In the following we consider particle configurations that belong either to C or to Cx [0,1]. 
This requires the obvious adaptation of all the definitions in Section 12.1.11 We perform 
the necessary adjustments in the following definition. 

Definition 5.2. Let {X,d) be a locally compact complete separable metric space. 

• A measure 6 on (X, Bx) is called a Radon measure if 0(B) < +cx) for any B E B\. 
It is called a particle configuration on X if in fact 9(B) E Nq for any B E B\. 

• The space of particle configurations on X is denoted by A/'*(X)j3 

• M*(X) is endowed with a measurable space structure by considering the a-algebra 
generated by the counting events on X. 


^Notice that in Section [2.1. 1 1 we asked the configurations in M(S x G) to be of locally finite allocation. 
Hence the extra index * in M* (X). This restriction was imposed in Section (2. 1.1 1 for simplicity of notation 
in the statement of Theorem 13.41 but it is unnecessary here. 
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Given a Radon measure i? on X we define the Poisson distribution as the unique 
probability measure on N'*{X) which satisfies 

7r^{{e e AA*(X) : e{B,) = ki ioT alU = 1,... ,n}) = Y[— ^ ^ 

i=l 




for all ki,..., kn E Nq, disjoint Bi,..., Bn E and n G N. 

• A random particle conhguration on X is called a Poisson process with intensity d 
if it is distributed according to 


In the following, we will often identify a given random cylinder conhguration V with a 
birth-and-death process on S' x G through its time sections; if for each t G M we dehne 
the random particle conhguration Vt E ^/{S x G) by the formula 

H({ 7 a:}) := #{C G V : basis{C) = 7 ^, and he <t <hc + Ic} 

for every 'jx E S x G, then V = {Vt)t&R constitutes a birth-and-death process on S' x G. 
From this point of view, we interpret any cylinder (y^., t, 1) as an animal 7 born at time t 
on location x with a lifetime of length 1. 

Now, consider a Poisson process IT on C with intensity measure (j),^ ■.= u x e~^^C x 
where C is the Lebesgue measure on M and is the exponential distribution of mean 1. 
We call n the free proeess, since it is a non-interacting birth-and-death process on S' x G. 
It is stationary and has ^ as its invariant measure. The process /C will be obtained as 
an appropriate thinning of 11. However, to properly conduct such thinning we need to add 
an additional component to H; to each cylinder in H we attach an independent uniform 
random variable, which we call its flag. Each of these hags will be used to determine 
the success of its cylinder’s attempted birth in the dynamics. More formally, we dehne 
the flagged free process H as the Poisson process on C x [0,1] with intensity measure 
01/ := 01/ X -^[ 0 , 1 ]- Foi’ any {■jx,tfl) G H we denote its corresponding hag by F^'jxytfl). 
Thus, elements of H are simply pairs of the form (G, F{G)) with G E C. Finally, we dehne 
the thinned process /C by the formula 

(5.2) K = {(7x^,0 e n ; F{'yxfl,l) < M{'yx\Xt-)} 

where, for G S' x G and E Af{S x G) we use the notation M{'yx\f,) ■= 

Observe that the self-referential nature of the thinning rule in fl5.2p could keep the process 
K, from being well-dehned. Indeed, let us introduce some dehnitions that will help us give 
further details on this matter. 


Definition 5.3. 

• Given G, G G C we say that G is a first generation ancestor of G and write G ^ G 
whenever 

hasisiC) hasis{G) and < be < 

where ^ is the impact relation dehned in (12.31) . We will denote by V{G) the set of 
all hrst generation ancestors of a given G G C. 

• For G G C we dehne AflG) := the restriction of H to V{G), and for n G N 

we set 

A+i(G):= IJ Mi(G). 

C&A„{C) 

We dehne the clan of ancestors of G in H as 

A{G) := IJ An{G). 

nSN 
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• For t G M and A E Bg let us define the clan of ancestors of Ax G at time t as 

A\A xG):= IJ A*JA X G) 

neNo 

where ^o(A x G) := {G G 11 : basis{G) eAxG, be ^'t< be + Ic} and for n G M 

X G) := IJ An{G). 
e&Al{hxe) 


Let us return to the discussion of the well-dehniteness of the process /C. Notice that if 
we wish to determine whether a given cylinder G = {'jx, t,l) eH belongs to /C or not then 
hrst we need to specify the conhguration /C^- in order to evaluate whether the condition 
in fl5.2p is satished. To be more accurate, due to Assumptions 12.151 we will only need to 
specify JCt- inside the set Hence, recalling Dehnition 15.31 we see that to determine 

the fate of G we must first determine the fate of, in principle, all of its ancestors in Ai{G). 
But this itself involves determining the fate of a second generation of ancestors, A 2 (G). 
In general, to determine if G belongs to /C one may need to study the fate of every cylinder 
in the clan of ancestors of G. If A{G) were to span over an inhnite number of generations 
then it may be impossible to decide whether to keep G or not and, hence, the process K, 
would not be well-dehned in this situation. On the other hand, if for every cylinder G G H 
the clan A{G) spans only over a hnite number of generations then /C would effectively be 
well-dehned. Indeed, since M{^x\A^1-) = Af(7a,|0) for any {'yx,tJ) ^ H with no ancestors 
preceding it, the fate of every cylinder in the last generation of ancestors of a given G G H 
can be decided upon inspecting their respective hags (and nothing else). By proceeding 
one generation at a time, the fate of all their descendants, including G, can be determined. 
More precisely, take G G H and dehne 

Nc = max{n G N : An{G) ^ 0}. 

If Nc < +CX 0 set 

■.= e Ano(C) : < M(%m 

and for 1 < i < Nq — 1 inductively dehne 


/G(G) = /0+i(G) U {(%,r,0 G Ai{G) : F(%,r,0 < M(%|iF,+i(G))}. 


Then the cylinder G G H is kept in K, if and only if 


F(G) < M(%|iFi(G)). 


This algorithm decides whether G is kept in hnitely many steps. Therefore, we have that 
for /C to be well-dehned it suffices to have Nc < +cxo for every G G H. As it turns out, 
the heavy diluteness condition of Dehnition 13.11 is enough to ensure this. 


Proposition 5.4 ( |211 Proposition 9.3]). If (p, Lf) is heavily diluted then almost surely 
the clans of ancestors A^{A x G) are hnite for every f G M and A E Bg. 

It follows from the discussion above and a straightforward computation that if {u, H) 
is heavily diluted then /C is well-dehned and its inhnitesinial generator is given by fl5.ip . 
Furthermore, Proposition 15.41 implies the following facts: 

i. The dynamics dehning /C loses memory of the initial condition, a fact which implies 
uniqueness of the invariant measure. Furthermore, since /C is a stationary process 
(due to its time-translational invariant construction and the stationarity of H), 
this implies that for each f G M the conhguration Kt is distributed according to 
the unique invariant measure of the dynamics. 
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ii. /Cf/Ci as A S' for all t G M and rj G Af{S x G) of finite if-interact ion range, 
where JC^ is the process obtained by rnnning the dynamics in the volume A with 
boundary condition rj (and with AE^ replaced by in the dehnition of M, 

see I5.2p . Each process is again stationary and a straightforward calculation 
with its generator shows that /xai?? is its unique invariant measure. From this and 
the local convergence stated above, we conclude that the invariant measure of /C 
is the unique gas measure of (z/, H). 


6. Proofs of Theorem 13.41 and Corollary 13.51 

In this hnal section we give the proof of our main result. Theorem 13.41 Following the 
proof, we discuss some relaxations of its hypotheses and then conclude the section by 
showing how Corollary 13.51 follows from Theorem 13.41 

6.1. Proof of Theorem 13.41 We divide the proof in three steps. 


6.1.1. Uniqueness of the gas measure for each model (n^, H^). Let us begin by considering 
a Poisson process If on C x [0,1] with intensity measure := u x e~^^C x £r+ x W[ 0 , 1 ]. 
Notice that, by the nature of the intensity, all cylinders in If have multiplicity one. Thus, 
it makes sense to dehne for each £ > 0 the e-discretized process If (or simply e-process) 
by the formula 

(6.1) := {(7^,t,s,M) G C X [0,1] : (7„,f,s,M) G H}. 

Observe that 11^ is a Poisson process on C x [0,1] with intensity x C x £r+ x U[0,1]. 
Moreover, fl 6 .ip establishes a one-to-one correspondence between cylinders of fl and fl^. 
Thus, in the following we write C^- to denote the e-cylinder in fl^ which corresponds to 
the cylinder C G If, i.e. if C = (y^,, t, s) then we set = ( 7 ^, t, s). 

Now, by the proof of Theorem 13.21 we see that it suffices to show that for each e > 0 
^o,h®(a X G), the clan of ancestors at time 0 with respect to the Hamiltonian prescription 
and with underlying free process H^, is hnite almost surely for all A G To see this, 
notice that if we dehne an impact relation by the condition 


% lx ly e V (7^) 

and for A G we consider the corresponding clan of ancestors M°’''"(A x G) given by 
then for each e > 0 we have 

( 6 . 2 ) A°’^"(AxG) CE,(A°’^(AexG)) 

where A^ = {x G S' : d(x, A) < e} and for P C C we set 


Il,(r) = {(7^,t,s)GC:(7.,t,s)Gr}. 

By the proof of Proposition 15.41 (see [21] for details) we see that ’^ < 1 implies that 
almost surely M°’'^(A x G) is hnite for every A G Hg, so that from fl6.2p we conclude that 
A^'^^{A X G) must be hnite as well. This shows that {A, H^) has a unique gas measure. 
Notice that our argument does not imply that the model {A, H'^) is heavily diluted itself. 
Indeed, the only available estimate 


diYx) 




q{ly)dA{Ay) < 


-,-AE 


Q{Tx) 


'U(7T 


Qily)dl^\ly) 


for every G S' x G is, in principle, not enough to show that < 1 implies 


< 1 . 
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6 .1.2. Construction of the coupling We need first the following technical lemma. 

Lemma 6.1. Snppose that N C J\f[S x G) is a a-local vr'^-null event. Then 

(6.3) P(n.t £ N for some f G M) = 0. 

Proof. If iV = IJneN fo^ ^ seqnence (A^n)nGN of local events, then it snffices to show 
that for each n G N one has 

Pijit £ Nn for some f G M) = 0. 

Thns, let ns hx n G N and consider A E Bg such that iV„ G Paxg- Lef us observe that 
since is J^Axc-rneasurable we have 

(6.4) IIj G Nn G Nn- 

Now, let Ma be the set of cylinder conhgurations 9 G x G x M x M+) such that: 

i. 9t{A X G) < +CXD for all f G M. 

ii. Every cylinder in 6 has a strictly positive lifespan, i.e Ic > 0 for all C E 9. 

Notice that if 0 G Ma then the process 9 = {9t)t&. is piecewise constant. In particular, 
for 9 E Ma we have 

9t E Nn for some f G M 9r E for some r E Q. 

Now, the choice of the intensity measure (fi, plus the fact that i^(A x G) < +cxo yields that 

-P(nAxGxRxR+ £ Ma) = 1. 

Hence, from fl6.4p we obtain that 

F(nt G Nn for some t G M) = P((nj.)AxG ^ fo^ some r G Q) 


< P{Ur ENn) = J2 = 0- 

rSQ rSQ 


□ 


Now, we construct the coupling (.E^)e>o by considering the processes 
/C^ = {(7|,t,s) G n ; F(7^,t,s) < M^{YJ^.)} 
where for each E S x G and ^ G J\f{S x G) we dehne 

A/r(7,|0 := 

By the arguments in Section [SI each process /C^ is stationary with invariant measure 
Thus, for each e > 0 we may dehne Z'^ := /Cg. We need to check that for any B E B^^q 
there exists Eb > 0 such that for every e < Eb 

(6.5) Z% = D, (4,-(B)) ■ 

To see this, we notice that if N denotes the dynamically negligible event in the statement 
of Theorem 13.41 then by Lemma [6.11 and since /C° is a thinning of H we can assume that, 
with the exception of a set O of realizations of H with zero probability, for every G G H 
we have 

( 6 - 6 ) i^bc) I{basis{C)) ^basis(C) ^ N . 


What we will show in fact is that, for almost every realization of the process H outside of 
this exceptional set O, for any given B E Bg^Q if £ is sufficiently small then (16.51) holds. 
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The heart of the proof is contained in the next lemma, which states that clans of ancestors 
stabilize for e small. 


Lemma 6.2. Fix A E and let O be the exceptional set above. Then in there exists 
almost surely (a random) £o > 0 such that for 0 < e < £o 

(6-7) ^De(AOT(AxG)) ~ (^!40 T(AxG)) ’ 

i.e. any C E x G) is kept in IC^ if and only if Cg is kept in /C^ for every 0 < e < £o- 

Proof. Dehne K := max {n G N : A^^{A x G) ^ 0}, which is hnite since < 1. 

Notice that for every cylinder G E {A x G) and e > 0 one has 

G,ElC^ ^ F{G) < M%basis{GeM. 


Thus, since 0 + (5basis(G) ^ by fl6.6l) and I{basis{G)) = 0 by dehnition of K, by recalling 
that 



for every 'fx ^ S xG and rj E A/{S x G) with G N'^, we deduce that for e (randomly) 

small enough 

^DeiA^iAxG)) ^ (^A^''(AxG)) ’ 

except perhaps if F{G) = M^{basis{G)\0) for some cylinder G E A^j^{AxG), a fact which 
can only occur with zero probability. In a similar way one may proceed with the following 
generations by induction, using at each step (16.6j) and (I6.2jl . to arrive ultimately at (16.7j) . 

□ 


Now, £x i? G Bg^Q and take A E Bg sufficiently large so that i?d) C A x G, where i?d) 
is as in Dehnition 12 .2UI By Lemma [6.21 applied to this volume A and the inclusion 

D;\K) GAxG 

valid for every 0 < e < 1, we conclude that (16.5j) holds for e small enough. 


6.1.3. Concluding the weak convergence gd ^ gP. The weak limit now follows from fl6.5p 
by the dominated convergence theorem and the following lemma. 

Lemma 6.3. Let us suppose that (^^'^^)e>o ^ A/'(5' x G) satishes that for any B E Bg^Q 
there exists Sg > 0 such that ) for all e < eg- Then -A- vaguely. 

(-D) 


Proof. It suffices to show that for each compact set K C S x G and 5 > 0 there exists 
£o > 0 small enough such that G for all 0 < £ < £o- But if given B E B%^q 

we dehne 

PB ■= ^ min{d5xG {ixAy) ■ lx ^ 

and take 


^0 Pk^s) a £K(s) a 5 > 0 

where and are the sets from Dehnition l2.201 then G (^^°^)a ',5 for every e < Sq, 
since: 
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i. From the definition £K(s) and the fact that K C K{^s) we see that every G 

is of the form = D^(jx) for some jx £ Moreover, since s < p^is) there is 

at most one such so that, in particular, their multiplicities must be the same. 
Thus, the application p : —)■ [,^] given by 

pi'yx\i) = {D7\'yx^),i) 

is well-defined and injective. Since e < 6 we obtain that 

ii. By definition of Sr^s-) and the fact that K ^{K(^s)) we see that for every 

7a; G there exists G such that D^ipfx) = Furthermore, since 

£ < PkW there is at most one ya; G being mapped by to yi^^ so that, 

in particular, the multiplicity of yi*^^ in is at least equal to that of ya; in 

Therefore, the application p : —)■ given by 

p(yx,*) = (T>£(ya;),^) 

is well-defined and injective. Since e < 6 this shows that 

□ 

6.2. Some relaxations in the hypotheses of Theorem 13.41 Notice that for the proof 
of Theoreni l3.4l we only required the dynamically negligible event N to satisfy the following 
two properties: 

• is closed under the addition of particles. 

• P(Jlt G N for some f G M) = 0. 

Thus, if we call any event N verifying these properties a dynamically impossible event, 
then condition (ii) in the hypotheses of Theorem 13.41 may be relaxed by requiring that 
the event in question be only dynamically impossible. However, we should point out that 
in most cases of interest the original condition is already satisfied and it is simple to verify, 
so that there is no real gain from this relaxation. Furthermore, we also point out that in 
the statement of the theorem all infima and suprema may be taken to be essential in the 
measure-theoretical sense, arriving at the more general statement of Theorem 13.41 below. 

Theorem 6.4. Let be an approximation family and suppose that (z/^, is 

a family of diluted models such that: 

I. For every e > 0 the intensity measure P is given by 

;= I/O o T)-h 

II. There exist a z/°-null set M C S' x G and a measurable set N C J\f{S x G) closed 
under the addition of particles satisfying: 





" 

AE := inf 

inf 

inf ART(y|) 


£>0 

Aet3° 








ii. P{Ilt G N for some f G M) =0, where H is the Poisson process on C with 
intensity measure x e~^^C x where AE is as in (i). [Notice that AE 
(and hence H) depends, in principle, on the particular choice of N.] 

iii. For every y^, G M'^ and all p G J\f{S x G) with p + S-y^ G one has 

lim AE,^'(y^) = AEf{px)- 
£-)- 0 + 

iv. There exists for each y^, G S' x G a set V{px) ^ Psxg such that: 
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For every e > 0 one has the inclusion 

There exists a size function g : S' x G —)■ [1, +cxo) which verihes 


a 


_ 


. •= sup , 

7 xe 5 xG Iqnx) Jv(^^) 

Then the conclusions of Theorem 13.41 hold, namely: 

(a) Each model admits exactly one gas measure /i^ 

(b) As e —)■ O’*', we have the weak convergence 


,-AE 




< 1 . 




/• 


( 6 . 8 ) 


(c) There exists a coupling of the measures such that for any B G 

there exists (a random) > 0 verifying that for aW e < Eb 

Z% = D, (Zl-, 


■'B 


In particular, Z^ Z'^ with respect to the vague topology. 


6.3. Proof of Corollary 13.51 To conclude this last section, we show that Corollary 13.51 
is a particular case of Theorem 13.41 Indeed, notice that, since v, for every A E Bg 

we have that 7r)( with density given by 


(6,9) 


ZLla) = j,-(--(A>tC)->'(AxG)) 


dn 


A 


n 




dv 


a fact which can be deduced from the integration formula (see m Proposition 3.1]) 


r p-w(AxG) r / \ 

/ f{a)dnl = ^ ‘'"”( 4 . ■ ■ ■. 73 

Jj\f{AxG) • J{AxG)’^ \i=i / 


valid for all bounded J^axg- measurable functions / and measures t; on S' x G of S'-locally 
hnite allocation, where u" above denotes the n-fold product measure of v. In particular, 
the models (z/, H^) and {u, H^) are equivalent, i.e. they produce the same gas kernels 

where is the Hamiltonian prescription given by 

( diZ 



It is not difficult to check that for each e > 0 the pair (z/, H^) satishes Dehnition 12.151 
Moreover, since for every G S' x G it is possible to verify that /'^^({ya,}) = /^^({ya,}) 
and we also have that 

= Ae|; 

for every A G and y G M{S x G), we see that the family of models falls 

under the hypotheses of Theorem 13.41 for the approximation family "D = given by 

the identity operator on S' x G for each e > 0. The result then follows from Theorem 13.41 
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